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ABSTRACT 

Ammonia inversion lines are often used as probes of the physical conditions in the dense 
ISM. The excitation temperature between the first two para metastable (rotational) levels is 
an excellent probe of the gas kinetic temperature. However, the calibration of this ammo- 
nia thermometer depends on the accuracy of the collisional rates with H2. Here we present 
new collisional rates for ortho-NH3 and para-NH3 colliding with para-H2 {J - 0) and we 
investigate the effects of these new rates on the excitation of ammonia. Scattering calcula- 
tions employ a new, high accuracy, potential energy surface computed at the coupled-cluster 
CCSD(T) level with a basis set extrapolation procedure. Rates are obtained for all transitions 
involving ammonia levels with 7 < 3 and for kinetic temperatures in the range 5-100 K. We 
find that the calibration curve of the ammonia thermometer - which relates the observed exci- 
tation temperature between the first two para metastable levels to the gas kinetic temperature 
- does not change significantly when these new rates are used. Thus, the calibration of ammo- 
nia thermometer appears to be robust. Effects of the new rates on the excitation temperature 
of inversion and rotation-inversion transitions are also found to be small. 
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1 INTRODUCTION 

Since its discovery in the interstellar medium forty years ago by 
ICheung et al.] ( Il969h . ammonia has been widely used as a probe 
of the physical conditions in a variety of interstellar environments, 
ranging from p re-stellar cores, molecular clouds, to external galax- 
ies (see lHo ife Townes 1983, for a review). The peculiar structure 
of the molecule makes ammonia lines excellent tracers of the den- 
sity and temperature in these environments. NH3 is a symmetric 
top molecule, whose rotational levels can be denoted by two quan- 
tum numbers, the total angular momentum J, and its projection K 
along the C3,. molecular axis. Owing to the possible relative ori- 
entations of the hydrogen spins, two distinct species exists: ortho- 
NH3 (K = 3n, with n an integer; hereafter 0-NH3) and para-NH3 
(K + 3n; hereafter P-NH3). As both radiative and non-reactive col- 
lisional transitions do not change the spin orientations, transitions 
between 0-NH3 and P-NH3 are forbidden. Each of the rotational 
energy levels (with the exception of those with = 0) is further 
split into two sub-levels which can be denoted either by the inver- 
sion symmetry of the vibrational wave func tions or by the s ymme- 
try index e = ±1 (see Fig. 1 and Eq. 29 in lRist etalJll993ln This 
splitting is caused by the inversion motion of the molecule, and the 
corresponding inversion transitions fall in the range X ~ \ cm. Elec- 



* In this paper, we denote each level by its symmetry index e, with a + sign 
for e = +1 and a — sign for e = — 1. We refer the reader to the Fig. 1 of 



trie dipole transitions rules (A7 = 0, ±\,IS.K = 0) prevent radiative 
transitions between different K ladders to occur. Consequently, the 
lowest inversion doublets in each K iadder (i.e. with J = K) are 
metastable; they can be relaxed only through collisions. 

For this reason, the relative population of the first two 
metastable inversion doublets of p-NHs, Jk,i: = li,± and 22.±, de- 
pends solely on the kinetic temperature. Considering only the first 
three doublets, 22_+ and 2[ +, and assuming that the popu- 
lation of the 2i,± doublet is rn uch smaller than that of the 22.+ , 
IWalmslev & UngerechtsI h983l) showed that the excitation temper- 
ature between the two lowest doublets is given by the analytic for- 
mula: 
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iRist etaljfT993l) for an energy diagram of the molecule. 



where Tq is the energy difference between the first two metastable 
doublets (~ 41.7 K), T is the kinetic temperature, C(22 — > 2]) is 
the collisional excitation rate (averaged over the symmetry index e) 
between the = 22 and 2i rotational levels, and C (22 — » li) is the 
collisional de-excitation rate between the 22 and 1| levels. Thus, 
if one knows the T12 excitation temperature, one can derive the 
kinetic temperature of the gas, effectively using ammonia inversion 
lines as a "thermometer". 

Observationally, Ti 2 can be determined by observing the hy- 
perfine components of the li,_ — » li_+ and 22,+ — » 22- inver- 
sion transitions. The inversion doublets have indeed hyperfine com- 
ponents, which are due to the interaction between the electric 
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quadrupole moment of the N nucleus and the electric field gradi- 
ent created by the electrons. If one assumes that the excitation tem- 
perature of each hyperfine components (within a given rotational 
state) is the same, then one can derive the opacity t(1i) and t(22) 
of the 1 1 and 22 m ultiplets from the relative intensity of each hy- 
perfine component ( Barrett et al. 1977). T[2 is then calculated from 
the following formulae ( Ho et aL,1979i : ,Hotzel et al.,i2002 ): 



-T-o/ln 



9 T(22) 

20 T(li) 



(2) 



In order to "calibrate" the ammonia thermometer, i.e. to com- 
pute the kinetic temperature from the measured T12 excitation tem- 
perature, a good knowledge of the coUisional rates of NH3 colliding 
with H2 is necessary (see Eq.[T}. Although a large number of mea- 
surements have been made on the NH3-H2 sys tem, using in partic- 
ular double resonance (e.g. lDalv&Okj|l970h . crosse d beam (e.g. 
Schleipenetai]|l993l) and pressure broadening (e.g. Iwillev et al.l 



2002h experiments, laboratory data generally do not directly pro- 
vide state-to-state rate coefficients. As a result, radiative transfer 
models can exclusively rely on theoretical estimates. On the other 
hand, laboratory measurements are crucial to establish the predic- 
tive abilities of theory and, in particular, of the potential energy 

surfaces (PES). 

Followin g the pio neering work of iMorris et alj h973h . 

IWalm slev & Ungerechtj ([1983) performed statistical equilibrium 
calculations based on the theoreti cal NHs-He coUisional rates of 
lGreeij ( ll98(]h . lDanbv et alj ( Il988l) then determined more accurate 
coUisional rates for collisions of NH3 with p-H2(7 = 0) and used 
them to recalibrate the ammonia thermometer. The scattering calcu- 
lations of Danbv et al. ( 198 8) were based on the ab initio NH3-H2 
PES of lDanbv et alT ( fl986l) . This latter was subsequently improved 
to investigate propensity rules at selected coUisional energies for 
ammonia collid ing with both p- and 0-H2 dOffer & Flow^ll99d : 
iRist et all 1 1993 ). More recently, high accura cy ab initio calcu- 
lations have been performed by Mladenovic et al.l to ex- 
plore the topographical features of the NH3-H2 interaction. In the 
present paper, we present new coUisional rates based on the deter- 
mination of a new, highly accurate, NH3-H2 ab initio PES. Am- 
monia and hydrogen molecules are treated as rigid rotors and H2 is 
further constrained in the scattering calculatio ns to be in the (spher- 
ically symmetrical) para-y = state, as in iDanbv et al.l h988l) . 
Hence, the main difference between the present coUisional rates 
and those of Danby et al. ( 1988) arises from the PES. These new 
rates are then used to estimate the robustness of the ammonia ther- 
mometer and the excitation of ammonia lines under the condition 
that prevails in cold molecular clouds, prestellar cores and proto- 
stars (T < 100 K). For this, we compute the excitation of both 
0-NH3 and P-NH3 using a non-LTE radiative transfer code. The pa- 
per is organized as follows. In §2 we present both the new NH3-H2 
PES and the scattering calculations. Our excitation computations 
and the comparison with earlier computations are presented in §[3] 
and ^concludes this article. 



2 POTENTIAL ENERGY SURFACE AND COLLISIONAL 
RATES 

2.1 Potential Energy Surface 

The ammonia and hydrogen molecules were both assumed to be 
rigid. This assumption is adequate here because i) the investi- 
gated coUisional energies are below the firs t vibrational excitation 
threshold of ammonia (v2 = 932.4 cm" ' lRaiamakietai]|2004l) 



and ii) the corresponding collision time scales are much faster 
than the inversion motion of ammonia. Monomer geometries were 
taken at their gro und-state average values, as recommended by 
iFaure et al. The average str ucture of NH3 was deri ved from 

the high-accuracy calculations of iRaiamaki et al.l ( |2004|) : ^nh = 
I.95I2Bohrs, and HNH = 107.38°. The ground-state average ge- 
ometr y of H2 was taken a s thh = 1.4488 Boh rs (e.g. iFaure et al] 
l2005h . The conventions o f lPhiUipset all d 19941) were employed in 
defining the NH3-H2 rigid-rotor (five-dimensional) coordinate sys- 
tem (one H atom lies in the (x, z) plane). 

The NH3-H2 PES was constructed using the following two 
step procedure: (i) a reference PES was computed from a large 
set (89,000 points) of CCSD(Til calculations using the Dun- 
ning's correlation consistent aug-cc-pVDZ basis set; (ii) this ref- 
erence surface was calibrated using a complete basis set (CBS) 
extrapolation procedure based on a smaller set (29,000 points) 
of CCSD(T)/aug-cc-pVTZ calculations. A CBS-type extrapolation 
was applied to the correlation part of the interaction energy and 
was performed using a two-point X^^ type extrapolation, where 
X is the cardinal number corresponding to the basis set, as de- 
scribed in ljankowski & Szalewicd ( 12005 ). The self-consistent field 
(SCF) contribution was not extrapolated but was taken at the aug- 
cc-pVTZ level. All basis sets were supplemented with midbond 
functions and al l calculations were counterpoise corrected as in 
Ijankowski & Szalewicz ( 2005,). The same strategv was recently ap- 
plied to H,CO-H, jTroscompt et al.ll2009l) . 

Grid points were chosen for 29 fixed intermolecular distances 
R (in the range 3 - 15 ao) via random sampling for the angular coor- 
dinates of H2 relative to NH3. At each intermolecular distance, the 
interaction energy was then least square fitted using a 120 terms ex- 
pansion for the angular coordinates, using Eq. (3) of IPhillips et al.l 
( 11994) adapted to the C3,. symmetry of NH3. This expansion in- 
cludes anisotropics up to = 11 and I2 = 4, where the integer 
indices /j and I2 refer to the tensor ranks of the angular dependence 
of the NH3 and H2 orientation, respectively. The CBS correction 
was fitted over a subset of only 46 angular terms with /[ < 7 and 
li < 4. We note that the expansion restricted to P-H2 (J = 0), in 
which all terms with I2 i= are eliminated, includes only 24 terms. 
The accuracy of the angular expansions was monitored using a self- 
consistent Monte Carlo error estimator. A cubic spline interpolation 
was finally employed over the whole intermolecular distance range 
and was smoothly connected with standard extrapolations to pro- 
vide continuous radial expansion coefficients suitable for scattering 
calculations. Tech nical details on the fitting strategy can be found in 
IValiron et al] ( I2OO8D . The accuracy of the final five-dimensional fit 
was found to be better than I cm"' in the long-range and minimum 
region of the interaction (R > 5 Bohrs). The accuracy of the above 
procedure was also checked against a moderate set (1,200 points) 
of "high-cost" CCSD(T)-R12 calculations which offer a direct way 
of reaching the basis set l imit value within a single c alculation, that 
is without extrapolation dNoga & Kutzelnigg||l994l) . The RMS er- 
ror between the final fit and the benchmark CCSD(T)-R12 values 
was found to be lower than 1 cm"' in the whole attractive part of the 
interaction {R > 6 Bohrs). We emphasize that the intrinsic accuracy 
of CCSD(T) calculations at the basis set limit is ~ 1 cm"'. 

Constraining H2 in its lowest para level (J = 0) is strictly 
equivalent to averaging the PES over the H2 rotational motion. 
The global minimum of this averaged PES lies at -85.7 cm"' for 



^ CCSD(T) stands for the coupled cluster method with noniterative evalu- 
ation of triple excitations. 
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Figure 1. Comparison bet ween the first four an gular expansion coefficients, 
v/,,„, , of tlie NH3-H2 PES lDanbv et alj il986l) with those from the present 
work, as a function of the intermolecular distance R. 

R=6.3 Bohrs, with H2 in an almost equatorial location, equidis- 
tant from the two closest H atoms of ammonia. A similar location 
was found for the global minimum of the NH3-He interaction, but 
with a significantly more shallow potential well at ~ -33 cm"' 
l lHodges & Wheatle\^l200ll) . The five-dimensional PES, including 
the anisotropy of H2, is of course qualitatively different: the global 
minimum, as deduced from our fit, lies at -267 cm"' for R = 
6.1 Bohrs, with H2 colinear with the C^y axis of ammonia at the 
nitrogen end. It is interesting to compare this result with the recent 
calculations of iMladenovic et aP ( l2008h : these authors found the 
global minimum of the NH3-H2 interaction at a similar location 
with a comparable, although significantly smaller, binding energy 
(-253 cm"' ). As their calculations were performed at a similar level 
of accuracy (CCSD(T) method and aug-cc-pVQZ basis sets), this 
difference most likely reflects monomer geometry effects. Detailed 
comparisons will be investigated in dedicated future works. 

Now, in order to compare th e present NH3-H2 PES with that 
employed by lOanbv et al.l jl988 l), we present in Fig. [T]a compar- 
ison of the angular expansion coefficient s v/.m, . The definitio n of 
these coefficients is given in Eq. (2.1) o f iDanbv et alJ ( Il986h and 
their values are listed in their Table 4j- Only the first four are plot- 
ted for clarity in the figure. Despite significant differences at short- 
range (R ^ 6 Bohrs), the overall agreement between the two sets 
of coefficients is quite reasonable. This comparison i) indicates the 
good quality of the ab initio calculations o flDanbv et alj l ll986h and 
ii) suggests moderate effects of the new PES on the dynamics, as 
shown below. 

Finally, it should be noted that the PES ofl Danbv et al] h986l) 
has be en previously checked^ against laboratory measurements: 
ISchleipen et al.l h993h and Iwillev et alj ( 12002 1 reported, respec- 
tively, symmetry-resolved state-to- state and broadening cross sec- 
tions. In both cases, a good overall agreement was obtained be- 

v;|„,| corresponds to in Eq. (2.1) o f lDanbvetal]il98d) . 



tween theory and experiment, suggesting the adequacy of the PES. 
Discrepancies were however noted, in particular strong propensity 
rules predicted by theory for NH3-P-H2 were observed in experi- 
ment to a much lower extent (see below). 

2.2 Scattering calculations 

The quantal theory for scattering of a symmetric top with a n atom 
or a s tructureless molecule like H2 (7=0) can be found in lOreenl 
( Il976h . The extension of the formalism to the scat tering of a sym- 
metric top with a linear mo lecule can be found in lOff^er & Flowen 
( I1990I) and iRist et al" I ( fl993l) . In the present work, calculations 
were performed using th e (nonreactive) scattering code MOLSCAT 
( lHutson&GreeiJl994lFl in which the extension to allow for the ro- 
tational structure of H2 is not yet implemented. Hence, the present 
calculations were restricted to collisions between NH3 and p- 
H2 (y = 0). Extension to p-Il2(y = 0,2) and o-H2(y = 1) is under 
way and is further discussed below. 

All calculations were performed at the close-coupling (CC) 
level. Inversion doubling was neglected and the inversion-tunneling 
wavefunction was simply taken as a linear combin ation of two 
delta functions centered at the equilibrium position (Greenl ll976l : 
lDavis&Bogg ^ ri9781 We actually tested this approximation on 
the NH3-He system by taking the inve rsion coordinate explic- 
itly into account, as done previously by IPavis & Bo gg j ('l98l|). 
To this aim, we employe d the high quality NH3-He PES of 
'Hodges & Wheatlevl yOoT), which does include the inversion de- 
pendence of the interaction. The inversion motion was found to 
have a negligible effect (less than 10 % ) on the rigid-body interac- 
tion potential and on the cross sections jScifoni et al .l2007h . as was 
concluded by Davis & Boggs (1981) from a lower quality potential. 
We note that van der Sandenet alJ il99^ obtained a comparable 
result for the NH3-Ar interaction. A similar conclusion is there- 
fore expected for the NH3-H2 interaction, although the inversion 
dependence of this PES is yet not known. 

We adopted the rotational constants A = B = 9.944116 cm"' 
and C = 6.228522 cm"'. The reduced mass of the system is 
1.802289 a.m.u. As the ortho- and para-levels of ammonia do not 
interconvert in inelastic collisions, these were treated separately. 
The coupled-channel equations were integrated using the modi- 
fied log-derivative propagator of Manolopoulos ( 1986). The radial 
propagation used a stepsize parameter STEPS=IO except at total 
energies below 30 cm"' where STEPS was increased up to 300 to 
constrain the step length of the integrator below ~0. 1 Bohrs. Other 
propagation parameters were taken as the MOLSCAT default val- 
ues. Calculations were performed for collision energies between 
~ 0. 1 cm"' and 700 cm"'. The energy grid was adjusted to re- 
produce all the details of the resonances, with an energy step of 
0.2 cm"' up to total energies of 150 cm"' and 0.5 cm"' from 150 
to 300 cm"'. All calculations also included several energetically 
closed channels to ensure that cross sections were converged to 
within 7 - 8% for all transitions involving / < 3. Thus, at the high- 
est investigated energies, the basis set incorporated all target states 
with y < 1 1 and 12 for o- and P-NH3, respectively. 

Excitation cross sections for 0-NH3 are presented in Fig. |2] 
for rotation-inversion transitions out of the ground state of 0-NH3 . 
Prominent resonant features are observed in this plot. These are 
caused by both Feshbach and shape type resonances. It should be 
noted that only Feshbach type resonances are observed in NH3-He 

* |http : //mm . giss .nasa . gov/tools/molscat| 
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Figure 2. Cross sections for transitions out of tlie Oo,+ rotational level of 
0-NH3 into levels Jk ^ ( J ^ 3) as a function of collision energy. 
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Figure 3. Comparison betwee n rate coefficients f or collisional de-excitation 
of 0-NH3 by P-H2 (/ = 0) of Danbvetal.1 il986l) and those from tliis work, 
as a function of temperature. Only transitions towards the ground state Oo,+ 
are displayed. 



collisions (e.g. iMachin & Rouelj|2005h . This difference between 
He and H2(y = 0) reflects the deeper potential well of the NH3-H2 
PES, as discussed in Sect. 12.11 Resonances are found to signif- 
icantly increase the cross sections at low energy and, therefore, 
the rate coefficients at low temperature. For example, at 10 K, 
the NH3-H2 rate coefficient for the ground-state transition of o- 
N H3 lo,+ — » Oo,+ is a fact or of 20 larger than the NH3-He rate 
of iMachin & Roue3 ( l2005t) . We note, however, that this factor re- 
duce s to 2.5 for the tran sition 2o + — » O0.+ at the same tempera- 
ture. IWillev et all ( l2002h also reported significant differences (up 
to a factor of 4) between P-H2 and He broadening cross sections. 
It is also noticed that the cross section for the Oo,+ — » 33 + tran- 
sition is much lower than for Oo,+ — > 33,- . This propensity rule 
was already obs erved in earlier calculations iOffer & Flower! I990I : 
iRist et al]|l993h but, interestingly, it was found to be considerably 
weakened for collisions with 0-H2 (J = 1) and it was obse rved 
experimentally only to a slight extent jSchleipenetai]| 19931) . On 
the other hand, it was found to be preserved when including the 
7 = 2 state in the P-H2 basis set ( Off'er & Flower 1990; Rist et^ 
Il993h . This inclusion was also found to change the absolute values 
of the cross sections, at a few selected en ergies, by up to a factor of 
3 ( lOffer & Flow"erlll99d : lRist et al.lll993h . Its effect on the average 
cross sections and rate coefficients is however expected to be mod- 
erate, typically 20-30 %. This was indeed checked in the case of 
ND2H-H2 calculations employing the present PES (Scifoni et al., 
in preparation). As a result, the rate coefficients presented below 
are expected to be accurate within typically 30 %. 

Cross sections were integrated over Maxwell-Boltzmann dis- 
tributions of collisional velocities and collisional rate coefficients 
were obtained in the range 5-100 K for all transitions involving 
ammonia levels with 7 < 3 (the lowest levels with 7 = 4 lie at 177 
and 237 K above the ground states of p- and 0-NH3, respectively). 
Higher levels and temperatures were not investigated in the present 
work because collisional rates with 0-H2 are required in mod- 
els considering temperatures above 100 K. These collisional rate 



coefficients are made available in the BASECOlI and LAMB^i 
databases as well as at the CD^ 

In Fig. [3] downward rate coefficients are presented for 
rotation-inversion transitions towards the ground state. The reso- 
nant features are found to be completely washed out by the ther- 
mal average. The present results are compared with the data of 
Danbv et a l] l ll988h . As expected from the comparison of the ex- 
pansion angular coefficient s (see Fig.[T|l, the new rates agree within 
a factor of 2 with those of IPanbv et alj ( Il988l) . We see, however, 
that there is no particular trend, although the present rates are gen- 
erally larger 

In Fig.|4] we show the symmetry- averaged rates C(22 — > 2]) 
and C(22 — > 1|) that appears in Eq. (T}, as a function of the 
temperature. The new rates are found to be larger than those of 
lDanbvetal1 ( ll988h by typically 15% but to follow closely the same 
temperature dependence. As a result, minor modifications of the 
ammonia thermometer are expected, as shown below. 



3 NON-LTE EXCITATION COMPUTATIONS 

In order to estimate the effect of the new collisional rates on the 
calibration of the ammonia thermometer, we have computed the 
excitation of both 0-NH3 and P - NH3 u sing the large velocity gra- 
dient code of Ivan der Tak et al] (l2007h . We have used the 0-NH3 
and P-NH3 collisional rates with P-H2 pr esented in the previous 
section, as well as those from lPanbv et al.l ( 11988,) , fo r comparison. 
The l atter were taken from the LAMBA database dSchoier et al.l 
I2OO5I) . Collision between NH3 and He were neglected, because, 
in addition to H2 being more abundant than He by a factor of 5, 



^ http : //www . obspm . £ r/basecol/| 
* http: //www. strw. leidenurtiv.nl/~Tnoldata71 



http : //cdsweb . u- s trasbg . fr/| 



© 2009 RAS, MNRAS OOO.dHs] 



P-NH3 - p-H, (J = 0) 




20 40 60 80 100 

T.in (K) 

Figure 4. Comparison between the rate coefficien ts for collisional ex cita- 
tion and de-excitation of P-NH3 by P-H2 (/ = 0) o f lDanbvetal]il986h and 
those from this work, as a function of temperature. These rates have been 
averaged over the symmetry index e 
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Figure 5. Rotational temperature between the 1 1 and 22 metas t able l evels 
of p-NHj computed using the collisional rates of iDanbv" et al.l i 19881) and 
those from this paper, as a function of the kinetic temperature. The dashed 
curve show the rotational temperature computed from Eq.[T] 



NH3-H2 collisi on rates are typically a factor of 3 larger than the 
NHs-He rates jMachin & RoueffllioO^ . Energy levels, statistical 
weights and Einstein coeffici ents were taken from the JPL database 
for molecular spectroscopy IPickett et alj|l99§) . For the calcula- 
tions using Panby et al.. (1988) collisional rates, the first 24 lev- 
els of 0-NH3 and the first 17 levels of p-NHj were considered 
(corresponding to energy levels up to 416 and 297 cm"', respec- 
tively). For the calculations using the new rates, only the first 6 
levels of 0-NH3 and the first 10 levels of P-NH3 were considered 
(up to 1 18 and 1 15 cm"', respectively). In both cases, we have ne- 
glected the hyperfine structure of the molecule, i.e. we have consid- 
ered that each hyperfine level within a given inversion level corre- 
sponds to the same energy level. While this hypothesis will lead 
to an overestimate of the line opacity for optical d epths greater 
than a few, it is valid if the line is optically thin jPaniel et al.l 
We have therefore chosen a column density to velocity 
gradient ratio that is large enough for this approximation to be 
valid (jV(NH3)/(rfv /<jr) = 10""' cm~V(kms ~' pc"'); the same value 
adopted bv lWalmslev & Ungerech5 h983h ). 

Fig.[5]shows the excitation temperature between the li and 22 
metastable levels computed using both sets of collisional rates and 
as a function of the kinetic temperature. This excitation tempera- 
ture is obtained by summing the populations of the e = +1 and 
6 = -1 within each rotational state. A P-H2 density of 10'' cm"^ 
was assumed. On this figure, we also show the excitation tempera- 
ture computed from Eq. ([T}, i.e. assuming that only the first three 
rotational levels are populated. As seen on th is figure, the excita- 
tion temperature computed using the rates of Danbv et alj l ll988h 
and the one computed using the rates presented in this paper agree 
extremely well; both values differ by less than 2%. We also no- 
tice that for kinetic temperatures lower than 20 K, the excitation 
temperature is well approximated by Eq. (T](, but it underestimates 
it at larger temperature, because higher energy levels start to be- 
come populated. The good agreement between the kinetic temper- 



ature obtained using the rates o f lDanbvetal1 ( ll988h and those pre- 
sented in this paper can be simply understood by examination of 
Eq. (III. In this approximation, the excitation temperature depends 
on ln(l + C21/C23). Although the ra tes computed in thi s paper dif- 
fer by ~ 15% with respect to those of lPanbv et alj ( fl983) . their ratio 
(and a fortiori the logarithm of their ratio) differ much less. There- 
fore, the relation between the excitation temperature and the kinetic 
temperature - or in other words the calibration of the ammonia ther- 
mometer - appears to be robust. 

Fig. m shows excitation temperature of several P-NH3 inver- 
sion transitions, as a function of the density, for a kinetic temper- 
ature of 10 K and the same column density and line velocity than 
in Fig. [5] Spectroscopic data and critical densities are given in Ta- 
ble [T] For a multi-level system, the critical density can be defined 
(in the optically thin case) as the density at which the sum of the 
collisional de-excitation rates out of a given level is equal to the 
sum of the spontaneous radiative de-excitation rates: 

where the summation is done over the 7^,^, levels (with energies 
smaller than that of the Jf^j: considered). With this definition, the 
critical density refer to a level, and not to a transition. 

For densities lower than 10' cm"\ the excitation temperature 
of the Ii,_ — > I1+ and 22_+ — > 22_ inversi on transitions computed 
using the rates presented here and those of lPanbv et alJ show little 
differences. For these densities, collisional de-excitation is negligi- 
ble, and the excitation temperature of these lines are close to the 
background temperature (2.73 K). For densities much greater than 
the critical density (i.e. > 10* cm"'), collisional de-excitation dom- 
inates, and lines are essentially thermalized. At intermediate den- 
sities, the li,_ — > Ii + and 22,+ — > 22- line excitation temperatures 
predicted using the rates from t his work is sfi ghtly lower than the 
one predicted using those from IPanbv et al.l This is because the 
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Table 1. Frequencies, upper level energies and critical densities of the P-NH3 and 0-NH3 inversion and rotation-inversion transitions considered in this paper 
Spectroscopic data are from the JPL catalog. Energies are relative to the fundamental rotation-inversion level of each specie, i.e. li + for P-NH3 and Oo,+ for 
0-NH3. The critical densities are given for a kinetic temperature of 10 K. They are computed using Eq.[3] and thus refer to the upper of level of the transition. 
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Figure 6. Excitation temperature of the P-NH3 invers ion transitions com- 
puted using the coUisional rates of lPanbv et alj il988h and those from this 
paper, as a function of the H2 density. A kinetic temperature of 10 K is 
assumed. 



de-excitation rates from this w ork are smaller (by about a factor 
two) than those of 'Danbv et al. for these lines. As a consequence, 
the critical densities of the corresponding levels is greater than pre- 
viously estimated, and the transition thermalize at greater densities. 
Larger differences in the excitation temperatures of the 32 - — » 32,+ 
and 3i__ 3i + transitions - for which critical densities are a few 
10* cm"' - are seen. For example, at a density of 10' cm"', the 
excitation temp erature of the 32.- — » 32,+ transition computed us- 
ing the rates of iDanbv et al.l is almost a factor two larger than the 
one computed with the rates obtained here. From the observer point 
of view, this has no consequences because the energy of the upper 
level of the transition is 123.6 K. At low densities, for the kinetic 
temperature considered here, the fractional population of this level 
is extremely small, and the predicted antenna temperature is essen- 
tially zero. 

Fig. [7] shows the excitation temperature of the 0-NH3 33._ — » 
33,+ inversion transition as a function of the density, for the same 
column density to velocity gradient ratio than in Fig.|6l but - since 
the upper level of the transition lies at ~ 124 K above the ground 



10'' 10° 
n(H3) (cm-=) 



10' 



Figure 7. Same as in Fig.|6]for the 0-NH3 33 - - 33,+ inversion Une and for 
a kinetic temperature of 50 K. The dashed Unes indicate negative excitation 
temperatures. 

level of 0-NH3 - a kinetic temperature of 50 K The behavior 
of the excitation temperature is similar to that of P-NH3 inver- 
sion lines; at low density, it is close to the background temper- 
ature, while it is thermalized at densities greater than 1 0* crn ~'. 
We predict, in agreement with ^ Walmslev & UngerechtsI ( Il983b . a 
population inversion for densities ranging between ~ 4 x 10' and 
~ 6 X 10' cm"' (note that the range in which the inversion oc- 
curs is slightly different for the two collis ional rate sets). This pop- 
ulation inversion was first predicted by Walmslev & UngerechtsI 
(Il983h . and has been studied in detail bv lFlower et alJ (Il99(]l) . The 
corresponding maser transition has be en since observed in several 
star formin g regions, e.g. NGC 6334 1 ( Beuther et al.ll2007h . As ex- 
plained by Walmslev & UngerechtsI d 1983k the lower level of the 



* The computations with the rates from the present work do not include the 
levels above J = 3. To make sure that these levels can be neglected at a 
kinetic temperature of 50 K, we have computed the excitation temperature 
for the transitions shown on Fig. 7 using Dandy's rates, but without con- 
sidering the levels above 7 = 3. These were found to be quasi-identical to 
those computed when the levels above 7 = 3 are considered. 
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transition (33 +) is de-populated by collisions to excited levels in the 
K = Q ladder, while it is populated by radiative transitions from the 
33 - level. For densities greater than ~ 4 x 10^ cm"^, the collisional 
de-population rate is greater than the radiative population rate of 
the lower level, and the inversion occurs. For densities greater than 
~ 6 X 10^ cm"', excited levels in the A" = ladder start to become 
populated and populate the 33 + level collisionally. This limits the 
maser gain to moderate values; for a density of 1 x 10'' cm""*, we 

predict a negative opacity of only t 2. 

Fig.|7]also shows the excitation temperature of the Oo,+ 
rotation-inversion transition of 0-NH3. This line, at a frequency of 
~ 572.5 GHz, was firs t detected towards OMC-1 with the Kuiper 
Airborne Observatory jKeene et al.|[l983l) . It was also detected to- 
wards p-Oph A with Odin space telescope jLiseau et al .I2OO3I) , and 
it will be soon observable with the Heterodyne Instrument for the 
Far Infrared (HIFI) on board the Herschel Space Observatory. This 
line is found to thermalize at densities of ~ 10' cm"'. Once again, 
little differences between the excitation temperatures computed us- 
ing the rates of Danby et al. and those from this work are seen. No 
significant difference were found between the excitation tempera- 
tures of the other rotation-inversion transitions that will be observ- 
able with HIFI either. 



4 CONCLUSIONS 

We have presented new collisional excitation rates of P-NH3 and 
0-NH3 with p-H2(y = 0). With respect to older computations from 
IPanbv et al.l ( Il988b . the present rates were found to agree within a 
factor of 2. In order to investigate the effect of the new rates on the 
excitation of 0-NH3 and P-NH3, we have computed the excitation 
of these species under physical conditions that are typical of dense 
molecular clouds, prestellar cores, as well as the outer envelopes 
of embedded protostars, using an LVG code. We found that the ex- 
citation temperature between the 1 1 and 22 levels computed using 
the new rates is almost identical to that computed using older rates 
at the low temperatures considered here 50 K). Thus, the cali- 
bration of the ammonia thermometer appears robust. The effect of 
the new rates on the inversion transitions (at cm wavelenght) or the 
rotation-inversion transitions that will be observable with Herschel- 
HIFI are also found to be small. 

Future works include extension of the present calculations to 
P-H2 (7 = 0, 2) and 0-H2 i J = 1) as well as to higher ammo- 
nia leve ls and kinetic temp eratures. Compa risons with double res- 
onance jPalv & Okalll970h . crossed beam dSchleipen et al 
and pressure broadening dWillev et al.ll2002h experiments will be 
also investigated, with the objective to establish the predictive abil- 
ity of the present PES and to distinguish between the predictions of 
the available PES. 
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